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■ In this paper the Solomon network is simulated by means of ID and 

2D Ising model with additional — not only geometrical — neighbors. 
A "social phase transition" at a non-zero Curie-like temperature is 
observed, also in one dimension. The critical exponent describing the 
behavior of the magnetization in the vicinity of the transition is also 
^ \ evaluated. 
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1 Introduction 



In many sociological models, the behavior of each person is influenced by 
the neighbors, or influences the neighbors. However, the neighbors in the 
workplace are different from the neighbors at home. Thus the workplace 
neighborhood and the home neighborhood can be simulated by using two 
lattices of A" sites each, the home lattice and the workplace lattice. In the 
workplace lattice, the people are numbered consecutively from i = 1 to i = N 
with helical boundary conditions. The same people also show up on the home 
lattice, but an different order which is a random permutation of the order on 
the workplace lattice. Thus each person has exactly one place on the home 
lattice, and each site on the home lattice is occupied by exactly one person, 
just as is the case for the workplace lattice. But the same person occupies 
two entirely different sites on the two lattices. 

Such a network of two lattices is called a Solomon network (SN) |TJ be- 
cause each person is shared equally by the two lattices, just as in the biblical 
story of King Solomon. 

Within each lattice we have the usual type of interaction (Ising model, n- 
vector model, Sznajd model, ...); added to it is the interaction of each person 
with its own image on the other lattice. Thus in a chain of people i with 
nearest neighbor connection, the variables at site i interacts with the variable 
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at sites i — 1 and sites i + 1 as well as with the neighbors of the site P(i) of 
the other lattice, where P is the permutation of the numbers i — 1, 2, . . . , N. 

A simplified version of this two-lattice model uses only one lattice, and 
defines the neighborhood of each site as being the nearest neighbors plus one 
randomly selected site anywhere in the lattice. For a chain the neighbors of 

1 are thus i — l,i + l,i + R where R is a random distance; these neighbors 
can be compared with those of the honeycomb lattice: % ± 1, i± L. Therefore 
it cannot be excluded that ordering is found also in one dimension. This 
simplified model is closer to small- world networks |2|, |3|. f|]. 

2 Model 

Here, we consider interaction between individuals in the same manner like 
interacting spins in the Ising model. The intensity of interaction among 
neighbors at both lattices is described by the coupling constant J. When J 
is positive, neighbors at sites i and j tend to have the same opinion or share 
the same interests (Si = Sj). The external magnetic filed H is then global 
trend or world-wide fashion and forces spins to have the same sign as H. 
Then the temperature T plays the role of "social temperature" and describes 
a degree of randomness in the behavior of the members of the community. 

Before the system evolution according to Metropolis algorithm is started, 
additional neighbors of the z-th spin are chosen. For a give spin, the neighbors 
are fixed during simulation. During the simulation for each site at the lattice 
(individual, member of society) the energy change AE = E new — E Q i d asso- 
ciated with a possible spin flip (mind change) from Si into —Si is evaluated, 
where the energy is given by the Ising form: 

E = -JY,S l S J (1) 

<i,j> 

The sum in Eq. (1) goes over all neighbors pairs. For a given temperature 
T (expressed in J / ks units) a spin is flipped if and only if the thermal prob- 
ability exp(— AE/ksT) is greater than a random number taken uniformly 
from the interval [0,1]. 

We consider ID and 2D Ising model with one (IN) or two (2N) additional 
randomly selected neighbors. The simulation is performed on N = 10 6 sites 
with helical boundary conditions and initially, at t — 0, all spins are pointing 
up (Si = +1 for % = 1,2, . . . , N). Later we calculate the magnetization 
m = Y,iLiSi/N. The time is measured in Monte Carlo steps (MCS) and 
one MCS is completed when all iV spins at the lattice where investigated 
(in type- writer order) if they should flip or not. The value of the "social 
magnetization" m is averaged over the last 10 4 MCS. 
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Figure 1: The temperature dependence of the absolute value of the average 
magnetization. Simulation takes 2-10 4 [MCS]. The value of the magnetization 
is averaged over the last 10 4 [MCS]. The systems contains one million sites. 



3 Results 

It is well known from Ising himself that for the ID case the Curie temperature 
Tq = 0. Introducing one additional neighbor somewhere at the lattice in ID 
case shifts the "social Curie point" towards T c ~ 1-9(1) (in units of J/ks) 
as presented in Fig. |I[ This rough estimation of Tc gives a little bit higher 
values than reported in H £§]. 

For a given temperature T the magnetization m(t) and its cumulative 
time average (m(r)) = Y7t=\ Hf=i S%/(Nt) are stable except of small fluctu- 
ations (as presented in Fig. ^). 

The original SN model of two interacting Ising chains gives higher Tq 
than the 1D2N case with two additional and random neighbors (see Fig. 
[[]). The difference between these two models is that in the latter case a 
spin may interact with itself (selfish/introvert individual) even twice (super- 
selfish/super- introvert individual) and each spin may interact with another 
one more than once (up to iV times in the worst case). 

For the ID IN case - - which corresponds to the P§kalski model with 
p = 1 § as well as to Hong-Kim- Choi model [[| - we have found the 
critical exponent (3 (describing the critical behavior of the magnetization in 
the vicinity of the transition) equal to 1/2 (see Fig. |3|) which is similar to 
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1D1N:N=10 4 , T=1.5[J/k B ] 
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Figure 2: Time evolution of magnetization m(t) and its cumulative time 
average (m(t)) for 1D1N case and T = 1.5 [J/ks]- The systems contains 
N = 10 4 sites. 



1 D1 N: L=1 , N iter =2 10 [MCS], T c =1 .9 [J/k B ] 




Figure 3: The magnetization (m) vs. (T c - T)/T c for 1D1N case. T c w 
1.9 [J/fcs]. The systems contains iV = 10 6 sites. The simulation takes 2 • 10 4 
[MCS]. The dashed line shows (m) oc ((T c -T)/T c ) l/2 . 
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the results of extensive Monte Carlo simulation presented in Ref. Q. We 
have observed that the average absolute value of the magnetization (\m\) for 
T ~ Tc decreases with the system size 10 < L < 10 4 , in contrast to Ref. [J3J. 
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